% Function handle om gebruikt te worden door ode45

function f = odefun(t,y)

global a;
global b;
global c;

f=zeros(3,1);
f(1)=y(2);
f(2)=y(3);
f(3)=a*(b^2+c^2)-(2*b+a)*y(3)-(b^2+c^2+2*a*b)*y(2)-a*(b^2+c^2)*y(1);
end